Distribution of Matter from Singularity with Spherical Symmetry Using Fick Diffusion 
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A model is presented allowing calculation of energy and matter distribution in the Universe after 
expansion from singularity without introduction of expansion energy. Beginning with Fick's law of 
diffusion, we solve the Bessel function for spherical systems during expansion, presuming isotropic 
matter distribution in the Hubble flow. This function can be substituted with the Associated 
Legendre differential equation and when solved we discover useful parameters for those of energy 
and matter densities, diffusion and temperature. Though we can follow the decline of matter density 
over spacetime, we cannot suggest a precise value at singularity. This model may be useful though, 
as a starting point for modeling star formation rates, SFR, galaxy formation horizon lengths and 
regions devoid of matter where tracing the matter density decline over time is important. 
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I. INTRODUCTION 

A prediction of the General Relativity Q , is the pres- 
ence of singularity in our past. It seems all energy, mat- 
ter and spacetime had once been confined within an in- 
credibly small point which is now a mathematical and 
physical problem. While high energy experiments are 
helping unravel many details of particle formation soon 
after singularity 2], cosmology might propose estimates 
for initial energy densities, particle and energy diffusion, 
temperatures, the nature of immature spacetime, Q, Q, 
Q Universe age, initial galaxy size and star formaton 
rates (SFR). 

A difficult problem of cosmology is the description of 
energy and particle distribution in spacetime after release 
of the Universe from singularity. Evidence from the cos- 
mic microwave background (CMB), observed from a time 
about 380,000 years after singularity, suggests that en- 
ergy distribution was extremely homogeneous 0] and re- 
cent estimates of the smoothness of this radiation present 
a Gaussian distribution of very small fluctuations Q- We 
also know that soon afterwards matter recombined into 
atomic H and He, quickly followed by the birth of dense 
stars, star clusters, galaxies and even black holes @. This 
early Universe, dense with matter, is thought to have 
been a violent era when huge stars devoured enormous 
amounts of matter ending "life" with incredible super- 
nova (SN) explosions, littering the surrounding space- 
time with many of the heavier elements we presently en- 
joy [9][nj]- This quick appearance of star-rich galaxies is 
puzzling and leads to suggestions that galaxy formation 
began before recombination jlllj . 

Much effort is being made in understanding the na- 
ture of galaxy and SFR by modeling matter distribution 
after recombination These semi-emperical models 
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are quite detailed and are usually based upon cold dark 
matter (CDM) and dark energy considerations (ACDM). 
These models indicate a maximum SFR between 11 and 
13 Gyr ago, very soon after recombination SFR 
models are first order dependent upon matter density 
decline and gas temperatures, along with several tru- 
ely empirical parameters. Unfortunately, SFR also seems 
dependent upon the nature of gathered data; dust con- 
tributes mightily to noise in the optical spectrum which 
complicates checking the models against relible data. 
Comparisons with radiotelescopic data have been made 
and the expected redshift dependence of SFR is observed 
but, significant observational technique dependence of 
data has also been reported ^4[- Enourmous dependence 
is also placed upon cold dark matter (CDM) which is gen- 
erally accepted but experimentally unconfirmed; perhaps 
existing as plentiful amounts of really exotic particles. 
CDM is thought to contribute much more to galaxtic 
gravity than shining matter, star remains, dwarfs and 
planets Dark energy may also play a role but is 

mathematically inconsistent and has changed from a re- 
pulsive to attractive force looking back to about z«5 ^| 
an important epoch for SFR. 

We have generalized the first moments of the Universe 
using some well-known laws to describe matter and en- 
ergy distributions back to near the beginning of time. 
We have used the general diffusion equation of Fick, as- 
suming a single point was divisible into a nearly infinite 
numbers of small and equally proportioned components, 
to solve early spacetime expansion. The isotropic dis- 
persion of matter allows assumpton of homogenous and 
isotropic, non-relativistic diffusion of identical and nearly 
innumerable components and we do not distinquish be- 
tween ordinary matter and CDM. We have solved this 
specific situation as radial spacetime and by using the 
Galilean metric for matter and coupled radiation, have 
been able to predict regions of high matter density and 
those devoid of matter. This model is useful for predic- 
tions of matter density decline and for solutions of galaxy 
formation, very old SN rates and SFR. 
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II. THEORY 

A. Diffusion Model of the Universe 

The general diffusion equation of Fick can be applied 
to the Universe as matter diffusion from Singularity to 
present 



v.Vp 



(1) 



where p is matter density and k is a diffusion constant 
and v is the expansion rate. 

This equation is well known in fluid mechanics as the 
advection-diffusion equations for diffusion processes. We 
apply the Galilean transformation to particle dynamics 
in the Universe, after the cessation of inflation; we sug- 
gest particle nights are best modeled as non-relativistic 
with radiant energy strongly coupled to matter both dif- 
fusing well below light speed. We will apply a Galilean 
transformation to simplify the non-relativistic expansion 
by rewriting diffusion in a moving coordinate frame. Our 
Galilean frame is a radial expansion with a radial velocity 
not necessarily identical to the Hubble flow 

r = r' — vt and t' = t 



and with the derivatives M 



V,=V r - = -+v.V, 



(2) 



Placing these into Eq. Q we get the density dependence 
upon the new time and radial distance scales 



d P 



KV£ip + V.Vr'P 



d P V72 



(3) 



For the diffusion of matter and spacetime from a spherical 
point into expanding spacetime let us suggest that p be 
a function of 3-dimensions and the variable temperature, 
which is itself a function of time, p = u(r' ,9' ,<f>')T(t'). 
This model for matter diffusion is the partial differential 
equation 

u— = Tk\I t ,u. 



Because the left and right hand sides are equal, both sides 
must also be equal to a constant 



1 dT 1 rv2 x2 

— — = -nVtu = -X 2 
Tdt> u r 



(5) 



which allows us the freedom of two independent equations 

(6) 



IF = 



-X 2 u. 



For subsequent presentation we will designate these equa- 
tions as "thermal relaxation" and "matter diffusion", re- 
spectively. We will now let t represent some function of 
time, so the solution of the "thermal relaxation" model 
is 



T(i') = T e 



-aV 



This model seems straight-forward as a grossly general 
description for times after cessation of inflation, but must 
be modified for the special conditions during that epoch. 

To begin the solution of the "matter diffusion" portion 
we will conform to the usual conventions where the V 2 
operator is the Laplacian and the V, del operator, is the 
total partial operator. Using these terms the expansion 
has the following forms in spherical coordinates 



i_ a 
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(7) 
(8) 



We will now allow r = v r , while assuming that 8 — 
and 4> = in Eq. That is, the averaged polar and 

azimuthal velocities of matter are zero when traveling in 
the new reference frame, consistent with the presumption 
of a homogeneous dispersion of particles beginning with 
t = 0. Such is also consistent with the very homogeneous 
nature of the CMB radiation. Since this operator in- 
volves angular and radial variables we shall separate each 
variable for a useful solution. The solutions of the angu- 
lar portions of the operators are the angular functions 
Yi, m (0', 0') presented in Appendix A and are well-known 
in quantum mechanics. To use the step function 



L 2 Yi 



l(l + l)Yi, m (6,< 



(9) 



We now rearrange to separate the variables 
1 dT 1 



(4) we suggest u(r', 9', cf>') = R(r')Y^ m (6' ', (j)') for the "matter 
diffusion" portion of equation 



Tdt' 



-kV2,u. 



kV 2 ,u 



-X 2 u. 



We will separate u into radial and angular dependencies, 
which in a compact for is 



[l] 
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We can now rewrite this equation after separation of 
variables and substituting with Eq. (jJJ as 



,2<?R{r') , o ,dR(r') 



dr 



(10) 

where k 2 — remembering that n is the diffusion con- 
stant. 

Let us now substitute for y — kr' and by placing y into 
Ea, ljl()|l we have the following equation 



2 d 2 R(y) | dR(y) 
dy 2 dy 



y 



+ (y 2 - 1(1 + 1)) R(y) = 0(11) 



We will now use a second substitution of R(y) = 
Z(y)/y/y. When we use this for R(y) in Eq. l(TT)l we 
are left with a solution which is in the form of the Bessel 
differential equation of the type 



y ,fm +/ _m +(v ,_ {l+l _ r)z(y) ^ (12) 



Two solutions of which are essentially spherical Bessel 
functions (some primitive solutions are presented in Ap- 
pendix B) which we now define as 



nJ l+1/ 0)^ ni(kr , )= [^0) 
/kr' V 2 Jkr' 



The physical solution for both of these is well-behaved; 
ji(kr'), so the solution of radial portion is the propor- 
tionality 



Ri(r') oc ji(kr') 
and the general solution of the Eq.© is now 



(13) 



;=0 m=-l 

Where we allow 4>' — (j), and 6' = 9, because assume the 
expansion from singularity behaves isotropically and has 
only the radial component as non-isotropic. The path to 
the final result below is presented in Appendix C, where 
the orthogonal solutions to the Bessel functions are pre- 
sented 



ji(k lm r')Yr(0, <t>)p(r', 9, 0; 0)dV (15) 



B. Substitution of the Hubble constant for A 

If one recalls r = r a(t) where r is the present Universe 
radius, and a(t) is the expansion parameter (Hubble's 
law) is 



dr 
~dt 



ro- 



da(t) 
~dF 



Br. 



(16) 



Dividing each side r, inserting the matter density p 
M/(4irr 3 /3) then taking the derivative we have 



dp 
~dt 

dp 
~dt 



3M dr 
(47rr 4 /3) eft 

1 dr 



-3p 



r dt' 



Then dividing by each side with p and inserting the Hub- 
ble Constant 



1 dp 
p dt 



— f- = -3H. 



One can derive the left hand side from our Eqs. [221 ■ 1241 
or El 



1 dp 
p dt 



-^1 = —A 



so that we arrive at a compact form for our constant 
connecting temperature decline with matter diffusion as 



A = 3H. 



(17) 



p(r',8,<P;t) = E A Un e- x l^ ]l{kln y )Yr{e ^ ) {lA) 



At some time in the past, the Galilean frame becomes 
unnecessary as small perturbations merge togther and 
the Universe is better described by the Hubble expansion. 



C. Path of light in the early Einstein Universe 

We suggest much of the energy of diffusion in the expand- 
ing Universe just after recombination can be considered 
uncoupled photons. These "excess" photons were sub- 
ject to the conditions of the Einstein Universe and can 
be treated with the usual metric 



^ 2 = -(i-^r^ 2 -rW 



2 ■ 2 

r sin 



z 2 dt 2 . 



The path of a photon, if collision free during the ini- 
tial moments of recombination, will be subject to space- 
time curvature in the usual manner with R a constant 
of spacetime curvature. In a Universe with truly homo- 
geneously and isotropically dispersed matter and energy 
we should find K to be quite large and the light path to 
be nearly, if not completely flat, even in the primitive 
Universe. We expect inhomogeneties in the polar and 
azimuthal dimensions were extremely small, contribut- 
ing only very slightly to the curvature of photon flight. 
Light in a primitive Universe without truly isotropically 
dispersed matter and energy may suffer a flight through 
curved spacetime and a slightly smaller value for R. Such 
radiation might be observed at later times as our CMB 
[l7) | with much of the tiny temperature ripples due to 
more recent interactions with matter and supernova ad- 
ditions but still containing information before recom- 
bination when more recent events are removed ItJ. 
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III. EXAMPLES OF CONDITIONAL 
SOLUTIONS 



The following examples might be considered models 
for our Universe approaching t — 0. These are based on 
a function describing a declining proportionality between 
matter density in spacetime and the time/radial distance 
from the origin. The reader should remember that A 2 = 
3H in the following examples which are based upon the 
proportionality 

r + vt 

The parameters H and k may be chosen from observa- 
tions and we should also remember that k 2 oc ~; the 
square of the matter density is inversely proportional to 
a diffusion constant. 



p{r',9,<f>;t) 



AC 



T(^p) sin(fcr) 
r(f) kr 



e- 3Ht Y °( 



3C sin(^)r(2-n) sin(fc/) c „ 3m ^ 



nk 



l-n 



kr' 



It seems that Eq. (|2T1) will be difficult to evaluate since 
the value of the pre-exponential portion is highly depen- 
dent upon the relative value of the radial distance as 
our lookback approaches singularity, while the exponen- 
tial portion, dependent upon Hubble expansion, is time 
dependent. One thing to note from the pre-exponential 
factor is that energy /matter density depends on the dif- 
fusion of both through spacetime, as reflected by k. This 
is also consistent with the properties of a gravitational 
field, for particle distribution in any gravitational field is 
dependent upon radial distances between objects. 



A. Spherical Symmetric Initial Matter Density 

With commencement of the Universe, matter distribu- 
tion was approximately isotropic and can be represented 
as p(r',0,4>;0) = p(r'). In this case Eq. can be 

reduced to only the radial integral with the help of rela- 
tionships in Appendix A, Eq. i|A3Jl 



A, 



4-k 2 



,„ o / r' p(r%(k lm r')dr' (18) 
o 



In this example we have choosen a total mass of incred- 
ible density, p(r',6, 0; 0) ~ l/(r' n ). By substituting into 
Eq. l|T8jl and attempting to calculate the initial matter 
density we find the result is highly dependent both on 
the radial distance and time approaching the Universe 
origin. The total for matter and energy is represented by 
C in the relationship below 



-4, 



I .in 



C^5 lfi 6 m .o / r' 2 —3i{k lm r')dr' 
v 71 " Jo r 



where y = ki m r' . This function may be solved using 
a suggestion from by allowing p = 3/2 — n and 

v = I + 1/2 for a solution of the integral for the radial 
distribution in Eq. Illot as 



t^J v (t)dt 



2T( 



2 



and with substiution for the integral becomes 

r(^+3) 



Ai, m = 4C 



Sl.oS m o- 



r(^) 



(19) 



(20) 



By substituting this into Eq. ll"14t and using S functions 
we may generate the general case for the initial asymp- 
totic appearance of energy and matter density in our Uni- 
verse 



1. Homogeneous and Constant Matter Density at t = and 

n = 

Let us presume that the simplest case for equation 
(1211 . n = 0, corresponds to a constant initial en- 
ergy/matter density p(r',9, 0; 0) = C/r' — C, and that 
energy /matter was homogeneously distributed at singu- 
larity, t = 0. We can now use Eq. ll2Tll by allowing 
p(r' , 8, 4>\ 0) = C =constant. 

p(r\0 M = 2C S ^^ S ^ k pe-^ 
irk kr 

Due to sin(0) = we are left with 



(22) 



This result seems both puzzling and logical; suggesting 
that matter density was not distributed at all outside 
singularity. This also suggests that homogeneously dis- 
tributed matter at singularity was impossible, at least 
with this model. Perhaps at the commencing of the Uni- 
verse, t — 0, constant matter density everywhere within 
Universe cannot be well described, except to suggest that 
matter density was everywhere else but at singular- 
ity, whatever that "else" was. Another interpretation, 
though contrary to common thought, might be that at 
the instant of singularity there was intially no energy or 
matter density, at all. This suggestion, though strange, 
may also be a valid estimate of singularity at origin. 
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2. The Matter Density Asymptote at t = and 
< n < 2 

We again borrow from 0], r(|) = an d &o,o = 
fc, Ao,o = \/3£T, and we shall evaluate our model at n = 1 
and n = 2 but not n = in Equation (12111 



we can evaluate once again 

? sin (rrm + §) T(l — 2m) sin(fcr') _ 3J?t 



p(r',9,<f>;t) = 2C- 



irk 



kr' 



= 2C , (-l)"T(l-2m) simV) e _ 3m 



7rfc 



p(r',6,<t>;t) 



c 2 S in(fcrO e _ 3gt; n=1 
c sm(fcrQ e _ 3H n = 2 



DO 



(25) 



(23) This means the energy/matter density, everywhere and 
at all times, is always infinite. This might be interesting 
but not believable. 



Here we have arrived at an estimate of the initial mat- 
ter density of the Universe approaching singularity with 
some real values for the pre-exponential portion of our 
model. As expected from logical considerations and the 
arguments of |2J| the total energy / matter is still incalcu- 
lable, except to suggest it approaches infinite. 



3. An example of the above with the boundary condition: 
n > 2, n —even integers 

Let us now presume for n > 2 = 2m, m = 2,3,4, 
by again substituting real energy /matter density into Eq. 

2m 



2TJ p{r',9,(j);0) = C/r' zm and 

sin {rrm) Y (2 - 2m) sin(fcr') _ 3Ht 
p(r , 9, <j>; t) = 2C r-r-^r ^^-e 



irk 1 - 2 ™ kr 1 
(-l) rn - 1 k 2m - 1 sin(fer') 



3HI 



T(2m - 1) kr' 



(24) 



Here, the very left-hand portion of the relationship is 
again finite and the exponential a function of t similar 
to the previous example of n = 1,2. The complete pre- 
exponential cannot be evaluated at r = though. The 
above solution may be useful as a model, for at large 
times, radial distances and with a realistic diffusion con- 
stant the matter density can be calculated. As time and 
radial distance approach the infinite, the energy /matter 
density declines to 0, one predicted fate of our Universe. 



4- An example of a singularity created with the boundary 
condition: n > 2, n — odd integer 



Let us presume for n — 2m + 1, m = 1, 2..., and with 
substitution of real matter density into Eq. ll2l)l that 



p(r',6,ct>;0) = C/r 



/2m+l 



If we develop this Universe 



along the manner of previous examples, spherical symme- 
try but with spacetime in the forward looking direction, 



B. Solutions of a Point-Like Universe at t = 

1. Squeezed in a little sphere of r = a 

We again examine the Universe at t = and allow the 
matter density to be a smooth function of the spherical 
spacetime density with p ~ 9 (a — r) , where 9 (a — r) is 
the Heaviside, step function. If a — r > it obtains the 
value of 1, otherwise it is 0, which truncates the r integral 
at r — a. So we will examine a sphere with radius of a, 
where the density, C, may be equal to M/(47ra 3 /3). Let 
us suggest, for this solution, that p(r', 9, 0; 0) = CQ(a — 
r') and substitute this into spherical symmetric case Eq. 
I|18H then evaluate over all space 

4k 2 r°° o 

A l>m = -^S ll0 8 m , / r' 2 p(r')j/(fc ;m r')dr' 



a i 2 /■ oc 

C— ^<M™,o / r' 2 e(a-r%(k lm r')dr' 
Jo 

-SlflSmfl / r 3i{h m r')dr' 



(47ra 3 /3) 



3Mk f sin(fca) acos(fca) 



^3/2 a 3 



k 2 



k 



Further progess may be made by evaluating the series 
expansion for small values of a which obviously demands 
large initial matter densities 



At 



(26) 

After substitution into Eq. (114B and evaluation of the 
S functions the result, in terms of the Hubble constant, 
is dependent on the familiar parameters of time, radial 
distance, matter and diffusion 



Mk 2 sin(fcr') . 



-3Ht 



(27) 



For the solution leading to Eq. il27Tl we used the first 
term of the expansion when we took the limit of Ai iTn 
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at a = for the shrinking sphere. The solution of Eq. 
lj2TI) is infinite initial matter density at singularity, in as 
compact and precise form as can be presently suggested. 
This solution is similar to our Eq. 12411 where the matter 
density approaches infinite as the diameter of the Uni- 
verse approaches 0. So we have obtained similar results 
from two differing initial conditions. 



2. Point-like in Radial Direction Using the Dirac Function 

Imagine the Universe at singularity being compressed 
to a point with a volume of zero. It is possible to repre- 
sent this mathematically by a special function; the Dirac 
delta function and we shall use this to represent a tiny 
particle by a single point. Let us assume at t = that 

p(r' ,8,4>;0) = C 2 — S(r'), where we do indeed have 
r' sin 9 

real energy /matter density at singularity. We have sub- 
stituted this condition into Eq. l|l"5)) and have evaluated 
as follows 



.4 



l, in 



9k 2 



c 



2k 



I in 



7T 



ji(k lm r')Y l m (0,cj ) )p(r\9,cj ) ;O)d 3 r' 



ji(ki m r')5(r')dr' 



o 

2-7T />7T 

< I I Y l m (9,<p)d9d<p 
'o Jo 



171,0- 



(28) 



The angular integral of the spherical harmonics can be 
evaluated as in Appendix A Eq. (|A4j) . Here, we may use 
real values for matter density, diffusion and the expansion 
rate of the Universe; ji(0) = (5;,o 21] and substitute into 
Eq. ifTHl and again using 5 functions 



(29) 



We now have another result of interest, arrived by taking 
a separate path, for at most any radial distance from sin- 
gularity, given the diffusion constant, time and expansion 
rate, we should be able to calculate the matter density 
of the Universe, in a seemingly straightforward manner, 
using Eqs. (1271 l20|i . 



IV. DISCUSSION 

We have presented a model for the distribution of mat- 
ter in the Universe beginning at singularity with de- 



pendence on time, Hubble flow, and initial matter den- 
sity. By choosing as our initial condition the coupling 
of energy /matter density with thermal relaxation and 
Galilean comoving coordinates, we have been able to 
solve spacetime expansion using Bessel's spherical func- 
tions. To do this we first presumed expansion in the ra- 
dial direction away from singularity but as spacetime ex- 
pansion progresses the radial dimension degenerates into 
the present situation of relative spacetime. By separating 
Eq. (01 into two equal portions we have also assumed a 
constant value for total energy and matter immediately 
after singularity. We have chosen to not introduce con- 
tinuous synthesis of energy after singularity, even with a 
term similar to that first presented by Einstein Q and of 
current interest [23; this tends to cloud straight-forward 
evaluation because some solutions display regions of dis- 
continuity [ifH . 

We have used our solutions to examine several models 
for the Universe close to singularity. Such tests are im- 
portant because models of this epoch must be inferred 
from radiation emitted at much later times, which has 
been modified by intervening epochs , and current high 
energy experiments. The examples from our model are 
consistent, to the first order, with current CMB data by 
being homogeneous and angular independent. 

In Sec. IIIII we attempted to calculate energy/matter 
density back to singularity, that is towards t = 0. We 
presumed the initial matter density as spherically sym- 
metric and used r~ n type functions with asymptotes at 
r = 0. This may be the simplest model one can use to 
describe singularity with respect to energy /matter den- 
sity. The general solution of these type of functions, Eq. 
l[2"4"l) may seem straightforward, but insertion of specific 
values for n, for instance, selecting either even or odd 
values for n > 2, leads to differing results depending 
on the pre-exponential factors, which in turn depend on 
the integer chosen for the solution. An example solution 
taking this approach is (l22l) . which is the special case of 
n = 0, appearing not radially homogeneous, a condition 
we think is required by current CMB data, but we sus- 
pect this means everything but singularity exhibited no 
energy /matter density - at least to an "outside" observer. 
Another example 1251 describes matter density as infinite 
everywhere at real values for the radius of the Universe, 
so this model is not useful. 

We have found, through example IIII A 21 the n = 1,2 
case, and Eqs. I12--SII that we may only estimate the matter 
density at singularity as incredibly large as the Universe 
dimensions approach 0. We have found that matter den- 
sity close to singularity is also very dependent on matter 
diffusion, as expressed with the parameter, n again in 
Eqs. (l23l) . Examining these equations at the extreme 
situation of no matter, our model is consistent with that 
of DeSitter, where spacetime is absolutely flat and infi- 
nite, for when k = oo all matter must disappear. So this 
model, like the DeSitter Universe, falls short of observa- 
tions. When another solution, Eq. I|25|). is evaluated, 
we find the Universe consisted of infinite matter density 



while displaying a non-zero radius. This seems impossi- 
ble to us because there seems no reason for the eventual 
breaking of this "zone of infinite matter density" and al- 
lowing, spacetime expansion, thermal decline and matter 
diffusion. So this model is interesting but not useful. 

Testing the examples in subsection IIII B II for a Uni- 
verse with diameter r > we found that as the Universe 
is compressed towards singularity we obtain a solution, 
Eq. 112 7JI , with some similarity to that derived previously, 
Eq. i!24fl , but with a pre-exponential term we may possi- 
bly trace as close to singularity as we wish. For derivation 
(1271) . energy and matter densities are highly dependent 
on the radius of the Universe, the initial matter density, 
spacetime expansion and diffusion but will be incredibly 
large and nearly impossible to calculate accurately at sin- 
gularity. When we use the Dirac S function, in section 
IIII B 21 we also derive a matter density, Eq. I|29J1 . both 
similar to Eq. 1(271 and highly dependent on the param- 
eters mentioned for that equation. These two solutions 
may allow us another avenue to estimate the initial en- 
ergy/matter density of our Universe if we can estimate 
values for the Hubble constant, matter diffusion and Uni- 
verse age with accuracy near singularity. 

We can present an interesting example of introducing 
a perturbation, which we shall denote as x, into Ea . l(29"ll . 
We shall substitute with x by using the proportionality 
introduced in Section III, as the Galilean metric. If we 
let x — kv/H we can recast equation l12i)t in terms of the 
Hubble constant, time, distance and inversely propor- 
tional to a general matter difussion constant, compactly 
written as 



P_ _ sin(fcr + X Ht) _ 3Ht 
Pa (kr + xHt) 



(30) 



This equation allows us the flexibility that some fraction 
of matter diffuses with a slightly differing velocity than 
the Hubble flow. When we force x to be zero, for the 
unperturbed case, we show in Figure 1 that we can follow 
the increase of the matter density with lookback time 
as a function of the Ht while selecting for convenient 
values for the mean diffusion contant. We cannot find 
evidence, from this plot, for an epoch of discontinuity 
or "jerk" as proposed by Reiss and coworkers [22]; such 
must come from models using an explicit cosmological 
constant This situation of Fig. 1, also predicts 

epochs of negative matter density but only for special 
situations returning towards singularity. So this model 
appears much too general to indicate zones of lacking 
matter or rich in matter density. 
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0.25 
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x=o 



kr=0 

kr=Jt/2 

kr=3;i/4 

kr=7t 

kr=37i/2 



' t ! + / ++ 4 ' + + ' + + ' t+ ' ++f++ ' 4 



0.5 1 
Ht 



1.5 



FIG. 1: The behavior of the common solution from present 
matter density traced back towards early matter density in 
Eq. (!-'->( )tl for selected values of kr at t — 



to the kr axis, regions devoid of matter are predicted. 
Likewise, other regions become enriched in matter with 
increasing Ht at the expense of neighboring spacetime. 



p/po 



X=5 




FIG. 2: The behaviour of the common solution from present 
matter density traced back to the early matter density (p/po) 
versus kr and Ht for X = 5 



We next display another case of this relationship as a 
surface in Figure 2, beginning at singularity with x = 5 
rather than 0. Here we see that matter density decays 
very rapidly with increasing Ht and that at significant 
times beyond singularity we see, by tracing regions close 



The model displayed in Fig. 2 may be useful for es- 
timates of the effects of perturbations within the primi- 
tive Universe, at non-relativitic matter velocities, effect- 
ing the horizon lengths which determine early galactic 
group volumes. Such perturbations, in more detailed 
forms, may be introduced into Eqs. i|27ll or l|29j) . for 
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instance in the manner of Eci, l|3()|) . Models including de- 
viations from a Universe of ideal symmetry and including 
perhaps several Galilean metrics, may be found useful to 
reflect the small inhomogeneities, observed today, across 
the cosmic microwave background. Such perturbations, 
which may have been extremely tiny just after singular- 
ity, may well have become the nuclei for gas clouds, the 
first stars, clusters and galaxies. 



APPENDIX A: EXAMPLES OF SOLUTIONS TO 
SPHERICAL HARMONICS 

The mixture of Associated Legendre and spherical 
functions can be presented in spherical harmonics as 



17 



2 ' + 1 M?(eo^)^ 
4tt (l + m)\ 



Spherical harmonics condenses all of the terms, through 
the relation between Legendre functions and the range of 

TO, TO = —I. ..I 

A few examples of the simplest solutions to spherical 
harmonics for l / j m (6', </>): 



1 1 

1 3 



2 V2n 



cos 9 



Y, 



±1/ 



T 2V^ SmWe 



= I./^(3cos 2 6i- 1) 

4 V 7T 



Y 



±1/ 



1 15 ■ n 

2 V 2vr 



±u 



: (Al) 

All of these solutions are related through conforming 
with the orthogonality relationship 
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)Y,V 



Sll'S-mm' (A2) 



27T />7T 



JO 



Y, 



l 



4vr (l + m)\J A<P 



x / P/™ (cos 9)d9 

JO 



21 + 1 r 

— 2ir5 m>0 / p(cos 

4tt Jo 



^2-7T(5 m ,o f ^ . di 



4tt 



r(z + i/2) 



Z! 



(A4) 



The result of the integral of the Legendre polynomial is 
consistent with the formula of 1 24] 



APPENDIX B: SOME SOLUTIONS TO THE 
SPHERICAL BESSEL FUNCTIONS 



Some of the primitive solutions to the spherical Bessel 
functions: 



3o(x) 
n (x) 

ni(x) 
h(x) 
n 2 (x) 



SIM 

X 

cos a; 
x 

sin x cos x 



x* x 
sin x cos x 



x 



,i ix . 3 

— ) sin x - cos x 

y x d x x 



, 3 1, 
v x 3 x' 



- sin x — y— ) cos x 



(Bl) 



27T />7T 



.70 



Y, 



i \ u , 1 



i sin 



/>27T />7r 



= 2^ [ [ 

JO JO 

xY t m (9, (/))Y Q (9, (/)) sin 9d9d4> 
— 2v / 7r5;.o^m,o (A3)This is based on a more general solution by Joachain (2 



r 2 ji{h m r)ji{ki m r)dr 



9k 2 



(B2) 
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APPENDIX C: CALCULATING THE 
EXPANSION COEFFICIENTS 

The problem is not quite solved for we must perform 
further evaluations to uncover the coefficients. This 
requires numerical solutions for each Ai_ m for combina- 
tions of l,m. To do this we must impose orthogonality 
conditions upon Spherical Bessel solutions and spherical 
harmonic solutions; more details of these solutions are 
presented in Appendixes A and B. The general way to 
perform this is to multiply these solutions by orthogonal 
polynomials containing different parameters and then 
integrate over space, using the volume integral 



Eq. JB2tl for the spherical Bessel solution and for 
spherical harmonics Eq. <|A2|) can be solved as 



1 = \ E E Ame- A ^ im ; W T2-£«<^ 
1=0 rn=-l lm 



Using Kronecker <5's to kill those orthogonal summations 
and continuing from I — ► /' and from m — + m' we have 



/>0O pin />7T 

..dV = / / / ...r' 2 sin 6d6d(j)dr' 
v Jo Jo Jo 



A sample is given below if we reconsider Eq. 11411 



oo I 



P (r',e,ct>;t) = E Ai.^^m^Yne,, 



1=0 m=-l 



to use orthogonality relations we multiply each side 
jl'm'iki'm'^Y™ (8, if) and integrate over space and then 
then denote left hand side with / of the equality 



/ = / ji'm'iki'm'r')^ (8,cf > )p(r',8,ct > ;t)d 3 r' 
Jv 

oo / r 
= E A l,me- X ''^ / jl'm'ih'm'r'W > 



1=0 m=—l 



y.]i(k lm r')Yr{8,4>)d\ 

OO I 

= E E a 

17™ (6», <f)Y{! v (8, <f) sin ddddcj) 



1=0 m=-l 



Jo 



• : : / r' ji{ki m r)ji(ki> m ,r')dr' 



x 



10 JO 



X 2 . .+ 7T 1 



9 k l 

Z "'I'm' 



These solutions are based upon the independence of time 
and now we introduce this with an exponential. So with 
the boundary condition of t = we find: 



Jo = / Ji'm'fa'm'r'W (8,4>)p(r',8, ( f>;0)d 3 r' 
Jv 



Al'm/ — 



TT 1 



(CI) 



Z K l>r„> 



This leaves us with the final result of 



2k 2 

Ai, ro = — ^ / j l {k lm r)Y l m (8,<t>)p(r',8,ct>;0)d 3 r'.(C2) 
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